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Abstract.  In  this  paper  we  show  that  there  exists  a  family  of  supersonic  neutral 
modes  for  a  compressible  mixing  layer  in  an  unbounded  domain.  These  modes  have  zero 
wavenumber  and  frequency  with  non-zero  phase  speed.  They  are  analogous  to  the  super¬ 
sonic  neutral  modes  of  the  compressible  vortex  sheet  found  by  Miles.  The  results 
presented  here  give  a  more  complete  picture  of  the  spectrum  of  the  disturbances  in  this 
flow. 
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1.  Introduction.  In  recent  studies  (Jackson  and  Grosch,  1989,  1990),  we  considered 
the  inviscid  spatial  stability  problem  for  the  compressible  mixing  layer  in  an  unbounded 
domain.  The  key  to  understanding  the  stability  characteristics  of  this  flow  is  the  under¬ 
standing  of  the  different  parameter  regions  for  which  various  types  of  instability  modes 
can  exist.  In  Figure  1  we  show  a  plot  of  the  sonic  phase  speeds  versus  the  Mach 
number,  where 


c 


c_=  Pt;  + 


VP, 

M 


Here,  M  is  the  Mach  number  of  the  fast  stream,  the  ratio  of  the  speed  of  the  slow 
stream  to  that  of  the  fast  stream,  and  P/-  the  ratio  of  die  temperature  of  the  slow  stream  to 
that  of  the  fast  stream.  These  curves  divide  the  phase  speed-Mach  number  plane  into  four 
regions.  If  a  disturbance  exists  with  a  Mach  number  and  phase  speed  in  region  1,  the  dis¬ 
turbance  is  subsonic  at  both  boundaries,  and  v.e  classify  it  as  a  subsonic  mode.  In  region 
3,  the  disturbance  is  supersonic  at  both  boundaries,  and  we  classify  it  as  supersonic- 
supersonic  mode.  In  region  2,  the  disturbance  is  subsonic  in  the  fast  stream  and  super¬ 
sonic  in  the  slow  stream,  and  we  classify  it  as  a  fast  mode.  Finally,  in  region  4,  the  dis¬ 
turbance  is  supersonic  in  the  fast  stream  and  subsonic  in  the  slow  stream,  and  we  classify 
it  as  a  slow  mode.  We  found  that  there  is  only  a  single  .subsonic  neutral  mode  for  two 
dimensional  waves.  Beyond  a  critical  Mach  number  the  Mach  number  at  which  the 
phase  speed  equals  that  of  a  sonic  wave,  the  subsonic  neutral  modes  are  transformed  into 
supersonic  neutral  modes  which  are  subsonic  in  one  ffeestream  and  supersonic  in  the 
other.  In  addition,  another  supersonic  neutral  mode  appears  at  M*  the  Mach 

number  at  which  the  sonic  speeds  of  the  fast  and  slow  streams  are  equal.  This  supersonic 
neutral  mode  has  the  opposite  behavior  than  the  previous.  That  is,  if  the  continuation  of 
the  subsonic  neutral  mode  is  supersonic  in  the  fast  stream  and  subsonic  in  the  slow 
stream,  this  new  mode  is  subsonic  in  the  fast  stream  and  supersonic  in  the  slow  stream. 

In  these  studies  we  did  not  find  any  neutral  or  unstable  modes  in  region  3  which 
were  supersonic  in  both  streams.  In  this  paper  we  now  show  that  there  exists  a  family  of 
neutral  modes  in  region  3  having  zero  wavenumber  and  frequency  with  non-zero  phase 
speed*.  The  approach  taken  here  follows  that  of  Papageorgiou  and  Smith  (1989)  and 
Papageorgiou  (1990),  who  considered  long  waves  in  the  wake  behind  a  flat  plate.  We 
show  that  these  modes  are  analogous  to  the  supersonic  neutral  modes  of  the  compressible 
vortex  sheet  found  by  Miles  (1958).  The  results  presented  here  give  a  more  complete  pic¬ 
ture  of  the  spectrum  of  the  disturbances  in  a  compressible  mixing  layer. 

2.  Formulation  and  Results.  We  consider  here  a  compressible  mixing  layer  in  an 
unbounded  domain  with  mean  velocity  U  (q)  and  temperature  T  (q),  where  q  is  the  simi¬ 
larity  variable  in  terms  of  the  downstream  distance  and  the  Howarth-Dorodnitzyn  variable. 
With  suitable  normalization,  U  varies  between  1  in  the  fast  stream  (q  =  <»),  and  <  1  in 


This  possibility  was  suggested  to  us  by  R.  Klein. 
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the  slow  stream  (Ti  =-oo).  The  temperature  T  varies  between  1  in  the  fast  stream  and 
Py  >  0  in  the  slow  stream.  The  results  presented  here  are  independent  of  the  detailed  form 
of  V  and  T . 

The  mean  flow  is  perturbed  by  introducing  two-dimensional  wave  disturbances.  The 
disturbance  equation  for  the  normal  velocity  perturbation  (j)  is  given  by  (Lees  and  Lin, 
1946) 


jU  -  c)V- 
G 


-c)(}>  =  0. 


0) 


where  primes  denote  differentiation  with  respect  to  the  similarity  variable  t],  and  G  is 
given  by 


G  =  TIT  -M^iU  -c)2j. 


(2) 


Here,  a  is  the  wavenumber  of  the  disturbance  in  the  downstream  direction  and  c  is  the 
corresponding  phase  speed. 

Since  we  are  interested  in  long  wavelength  perturbations  (a  — >  0),  the  appropriate 
expansions  are  found  to  be 


<{>  =  (t)o  +  a<t)i  +  a?-^2  +  •  •  • 

(3) 

c  =  Co  +  occj  +  a^C2  +  •  •  • 

(4) 

—  = - h  OC  G  1  -t-  ct^  G  2  +  ■  ■  ■ 

G  Go 

(5) 

where 


Go  =  T[T-M^iU  -Co)2], 


(6) 


Gi=-2ciM^Go^iU  -Co)T.  (7) 

Let  us  now  define  the  linear  operator  L : 


LQV)  =  iU  -  Co)H''  -  U'^  =  (G  -  Co)2 


-V  t 

G  -  Co 


(8) 


Substituting  the  above  expansions  into  (1)  gives  the  following  three  leading-order  prob¬ 
lems: 
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d 

dri 


(Go'Z.«t)o))  =  0. 


(9) 


-i 

0 


[L  ((t)i)  -  c  1  %])  +  ^  (C  j  L  (%))  =  0, 
dr] 


(10) 


(Go  *  [L  (<^)  -  C2  <|)o  -  c  1  <^l])  +  (G 1  fZ.  (4)1)  -  c  1 4>oI) 

+  ^(G2L(4)o))-(G  -co)4>o  =  0.  (11) 

dr\ 

The  appropriate  solution  to  (9)  that  is  bounded  for  large  Tj  is  given  by 

4)o  =  Ao(G -Co).  (12) 

where  A  o  Is  a  constant.  Based  on  this  solution  the  following  three  remarks  can  be  made: 

[1]  Since  solutions  in  region  2  of  Figure  1  are  subsonic  in  the  fast  stream  and 
supersonic  in  the  slow  stream,  we  see  that  4)0  is  a  solution  if  we  take  the 
leading-order  wave  speed  to  be 

Co=l.  (13) 


[2]  Since  solutions  in  region  4  of  Figure  1  are  subsonic  in  the  slow  stream  and 
supersonic  in  the  fast  stream,  we  see  that  4>o  Is  a  solution  if  we  take  the 
leading-order  wave  speed  to  be 

Co  =  Pt/.  (14) 

[3]  Since  solutions  in  region  3  of  Figure  1  are  supersonic  in  both  the  fast  and 
slow  streams,  we  see  that  4>o  Is  a  possible  solution  for  some  leading-order 
wave  speed  lying  in  the  range 

c_<Co<c+,  (15) 

where  c±  are  the  phase  speeds  of  the  sonic  neutral  modes.  Thus,  the  actual 
value  of  Cq  is  not  determined  at  this  order,  and  can  only  be  found  by  using 
the  higher  orders  terms  in  the  a  expansion. 

To  determine  the  leading-order  wave  speed  Cq  for  region  3,  we  consider  the  expansions  at 
the  higher  orders.  The  solutions  to  (10)  and  (11)  are  given  by 

4>i  =  — Cj  i4o  +  ^  1  (G  —  Cq)  f  T  ^r],  (16) 

Uu  -  co)^ 
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(1)2  =  -C2Ao+  CliU  ~  Co)j 


<I>1 


(U  -  CqY 


dr\ 


Ai(U  -  Cy) I 


(U  -  CoY 


dr]+  A2(U  -  Co)J 


(U  -  CqY 


dr] 


+  Ao(U j(U  -cofds]  dr\. 
(u  ~  Cor 


(n) 


We  note  here  that  both  (()i  and  <{>2  could  contain  another  term  that  is  proportional  to 
which  can  be  ignored  by  a  suitable  renormalization  of  A  q. 

Inspection  suggests  that  the  solution  diverges  as  Tj  ^  ±«>.  Thus,  these  expan^iio'ia 
will  break  down  when  1  "n  |  =  0(a“').  There  are  two  outer  regions,  one  in  the  fast  stream 
and  one  in  the  slow  stream,  which  will  be  considered  separately.  For  the  outer  region  in 
the  fast  stream,  the  appropriate  scalings  are 


Ti  =  a  '  z 


for  the  independent  variable,  and 


-  ^0  +  +  ■  ■  ■ 


(18) 


(19) 


for  the  dependent  variable.  Substituting  (18)  and  (19)  into  our  original  equation  (1)  yields, 
to  leading  order,  the  equation 


W  +  G$  00  =  0,  G$  =  m2(1  -  co)2  -  1. 


The  solution  for  outgoing  waves  is  found  to  be 


00  ~  ^0^ 


-i  z 


The  next  order  equation  is  given  by 

01  +  Gq  01  =  ^  00, 

whose  solution  for  outgoing  waves  is 


gi  =  -  Co), 


2ylG 


(20) 


(21) 


(22) 


(23) 


Expanding  the  solution  as  z  — >  0,  and  matching  with  the  inner  solution  (3)  as  q  — >  o®, 
yields  the  following  relationships 
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^0  —  ~  t^o)» 


B 1  —  —  c  1 A  0, 


(^0  ) 


For  the  outer  region  in  the  slow  stream,  the  appropriate  scalings  are 


+  ■  ■  ■  »  (28) 

and  z  is  defined  as  before.  Substituting  into  our  original  equation  (1)  yields,  to  leading 
order,  the  equation 

K  +  Go  ^0  =  0,  Go-  =  py  -  Cq?  -  liy].  (29) 

The  solution  for  outgoing  waves  is  found  to  be 


-i  z 


The  next  order  equation  is  given  by 


^1  +  Go  -  gi  ^0.  -  2ci A/^Py  (P(/ -  Co), 

whose  solution  for  outgoing  waves  is 

Expanding  the  solution  as  z  — >  0,  and  matching  with  the  inner  solution  (3)  as  q 
yields  the  following  relationships 


^0  “  ^o(Pt/  “  ^o)> 


Di  —  Cji4o, 


/l,=//lo(Pt/  -Co)^(Go)"‘'^. 
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2  ~  i 


(Go)''' 


p7’  c  I  (2  P/’  —  A/'  (p(/  —  c o)' )  ^  0- 


(36) 


The  two  conditions  (26)  and  (35)  on  A  i  yields  the  following  eigenvalue  equation  for 
Co'- 


pj-  [M2(Py  -  Co)'  -  Pc  ]  (1  -  Co)'*  Co)'  -  1  ]  (Py  -  Co)'*.  (37) 

This  equation  is  identical  to  (5.3a)  of  Miles  (1958)  if  we  re-express  his  result  in  our  nota¬ 
tion.  Miles  showed  (in  our  notation),  that  for  region  3: 

[1]  A  single  root  of  (37)  exists  for 


M  ^  Mt 


1  +  VPc 


(38) 


with  phase  speed 


Co  - 


Pc/  +  >/Pr 

1  +  VPc 


(39) 


which  we  classify  as  a  constant  speed  supersonic-supersonic  neutral  mode. 
Note  that  this  solution  is  independent  of  Mach  number,  and  corresponds  to  the 
phase  speed  at  which  the  two  sonic  speeds  c±  are  equal. 

[2J  A  double  root  first  appears  at 


A/  = 


(1  + 
r-Pt/ 


(40) 


with  phase  speed 


+  P/'' 

1  -h  p/'^  ■ 


[3]  There  are  three  distinct  real  roots  for 


M  > 


(1  +  P/'')''' 
1-Pt/ 


(42) 


One  of  the  roots  is  given  by  (39),  while  the  other  two  roots  must  be  found 
numerically.  For  the  special  case  of  P/-  =  1,  these  roots  are  given  by 


-7- 


1  +  Pt/ 


±  +  4  -  4Vm2(1-P(^)2+i 


1/2 


(43) 


The  root  which  corresponds  to  the  (+)  sign  we  classify  as  a  fast  supersonic- 
supersonic  neutral  mode,  while  that  which  corresponds  to  the  (-)  sign  we  clas¬ 
sify  as  a  slow  supersonic-supersonic  neutral  mode. 

In  Figure  2  we  display  the  real  roots  of  (37)  which  lie  in  region  3  for  pj;  =  0  and 
various  values  of  P^/-.  The  roots  are  the  neutral  phase  speeds.  For  completeness,  we  also 
display  the  neutral  phase  speeds  which  lie  in  the  other  regions  (Jackson  and  Grosch, 
1989).  The  classification  scheme  is  given  in  the  figure  caption. 

The  two  conditions  (27)  and  (36)  on  A 2  yield  the  result  that  Cj  =  0.  Thus,  we  have 
/l2  =  fli  =  Dj  =  0,  and  the  final  solution  in  region  3  is  now  given  by 


^  —  Aq{U  -  c 0) 


.  .  (1  -  Co)" 

1  -1- 1  a  — — 


I 


W  -  coY 


dr[  + 


f  -  CQ)^ds  1  dT[+  ■  ■  ■ 

^  iU  -  Co)2  ^  J 


(44) 


3.  Conclusions.  We  have  shown  that  the  compressible  mixing  layer  has  a  family  of 
supersonic  neutral  modes  in  region  3  of  Figure  !,  which  have  phase  speeds  which  are 
identical  to  tho.se  of  the  compressible  vortex  sheet  (Miles,  1958).  For  the  mixing  layer  the 
wavenumbers  of  these  modes  are  zero,  while  those  of  the  vortex  sheet  are  arbitrary.  One 
might  expect  this  since  the  vortex  sheet  has  no  natural  length  scale,  while  the  mixing  layer 
has  a  non-zero  thickness.  We  note  that  these  results  can  be  extended  to  three  dimensional 
disturbances  by  rescaling  the  Mach  number  by  cos0,  where  0  is  the  angle  of  propagation. 
Finally,  Artola  and  Majda  (1987,  1989)  showed  that  the  nonlinear  interaction  of  the  three 
supersonic  neutral  modes  of  the  compressible  vortex  sheet  can  lead  to  instabilities  of  the 
flow.  We  surmise  that  similar  behavior  might  be  exhibited  by  the  interaction  of  the 
supersonic-supersonic  neutral  modes  presented  in  this  paper.  We  are  currently  investigat¬ 
ing  this  numerically. 

Acknowledgements.  The  authors  are  indebted  to  R.  Klein  for  suggesting  that  a 
supersonic-supersonic  neutral  mode  might  exist  in  the  mixing  layer. 
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Figure  1.  Plots  of  the  sonic  speeds  c±  versus  Mach  number  showing  the  four  regions 
in  which  different  types  of  disturbances  can  exist. 


Figure  2.  Plot  of  the  neutral  phase  speeds  as  a  function  of  Mach  number  for  (i^  =  0 
and  (a)  (i-/-  =  2,  (b)  py  =  1,  and  (c)  py  =  0.5.  The  neutral  mode  classification  is;  (1) 
subsonic,  (2)  fast  (supersonic),  (3)  slow  (supersonic),  (4)  Cyv  =  1,  (5)  =  0,  (6)  con¬ 

stant  speed  supersonic- supersonic,  (7)  fast  supersonic-supersonic,  and  (8)  slow 
supersonic-supersonic.  The  sonic  curves  are  shown  as  dashed. 
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